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The Electric or Magnetic Polarisation of a Thin Cylinder of 
Finite Length by a Uniform Field of Force. 

By T. H. Havelock, M.A., D.Sc., Fellow of St. John’s College, Cambridge, 
Lecturer in Applied Mathematics in Armstrong College, Newcastle-on- 
Tyne. 

(Communicated by Professor J. Larmor, Sec. R.S. Received November 26,—Read 

December 13, 1906.) 

§ 1. Introduction. 

The problem to be considered is the calculation of the polarisation induced 
in a cylindrical rod of finite length by a uniform field of force parallel to its 
axis. 

The only formula previously given appears to be that of Green* for the 
normal component of the polarity at the surface of the rod. In the following 
investigation the method of solution of the general problem of induction by 
series is used, and expressions are given in §§ 4 and 5 for the mean moment 
of the rod and the induction across the central section. It is assumed that 
the susceptibility of the rod is independent of the applied force; consequently 
the formulae apply especially to a dielectric rod suspended in a uniform field 
of electric force, this being, in fact, the purpose for which the calculations 
were made. 

In the case ong thin rods in a magnetic field it has been assumed that 
the demagnetising factor might be taken as that of an ellipsoid inscribable 
in the rod; but this neglects the fact of the non-uniformity of the field 
within the rod, and various experimental investigations have shown that the 
simple hypothesis is untenable. In § 7 these experimental results are 
discussed in relation to the formulae obtained in the previous sections. 

§ 2. The Series Solution of the General Problem. 

Consider a body placed in a magnetic field whose potential is V, and let a 
constant tc denote the magnetic susceptibility of the substance. Let Q be 
the potential of the induced magnetism and <£ the total potential; then we 
have 

4> Q + V. 

* Green, ‘ Mathematical Papers,’ p. 106 ; also Maxwell, ‘ Electricity and Magnetism ’ 
vol. 2, § 439. 
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The components of the intensity of magnetisation I at any point are 
given by 


I = 



a 

07 /’ 





If dv be the inward-drawn normal at the surl..__ 
to be solved in finding Q may be expressed in the equation 

Q = J3(V±Q) ! 7 ?i 
J dv r 




( 1 ) 


where the integral extends over the surface of the body. 

The exact solution of this problem can only be found in a limited number 
of cases, as, for example, when the surface of the body is an ellipsoid. A 
method of approximating to a solution was given by Beer,* and has been 
further developed by Neumann,f liiecke,* and Wassmuth.§ The method may 
be regarded physically as one of successive superpositions in the manner 
of Murphy.|| Suppose, as a first approximation, that the potential V 
produces an induction Ii equal to — /c(d/dx, 3/3?/, 3/3 z)V ; then the distribu¬ 
tion Ii will have a potential Vi, which may in its turn be supposed to induce 
a distribution I 2 , and so on. The total induction will be given by the 
sum of the series I; and the success of the method depends upon the 
convergency of this series. 

Defining a series of functions by the relations 


-If?*! V,.®!!*,, 

J dv r J ov r 


SVj (Is 


and so on, 


( 2 ) 


the first form of series for the induced potential Q, given by Beer, was 

Q = 4tt*Vi + (47 tk) 2 Y 2 + (4™) 3 V 3 + ... . (3) 

But even if the intermediary potentials decrease in magnitude, as in most 
practical cases, this series for Q is convergent only for a limited range 
of values for tc. Other forms of series have been obtained which include 
Beer’s series, but give a wider range of convergency; two such series given 
by Kiecke and Neumann are equivalent and are expressed in terms of the 
functions Vi, V 2 ... by 

Q = pVi -r P 2 (Vi + Vo) + P»(V 1 + 2V 2 + V 3 ) + ... , (4) 

where p = 47t/c/( 1 -f 47 t/c) = 47r/c/ya. 


* Beer, 4 Elektricitiit mul Magnctisinus,’ p. 155. 

+ C. Neumann, 4 Das logarithm. Potential,’ p. 248. 
+ Riecke, 4 Wied. Ann.,’ vol. IS, p. 466, 1881. 

5^ Wassnmth, 4 Wied. Ann.,’ vol. 51, p. 367, 1894. 

|| Murphy, 4 Electricity,’ p. 93, 1833. 
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This series is certainly convergent for all practical cases in which the 
intermediary potentials form a decreasing series, for instance, in the case 
of an ellipsoid placed in a uniform field parallel to an axis ;* for then we 
find 

Y = -H* ; V! = iEn; V 2 = - ; .; 

where L is the usual constant associated with an ellipsoid and is less 
than 4 7r. 

Forming the functions for series (4) we obtain 

pVi = -LHa:; p 2 (V, + V,) = L) - -LHr; 

p 3 (Vi + 2V 2 +V 8 ) = ( ' 47r ~P V - LHa-, etc. 

Hence we have 

Q =«LHx (i -f \. 

ML M M~ J 

But this series is convergent for all values of *, and in fact gives the known 
exact solution 

Q 

We shall use, then, the series (4) in the following work, in which we consider 
the case of a long rod in a uniform field. 


§ 3. Series for a Finite Rod. 

We shall consider the case of a cylindrical rod of length 21, with flat 
circular ends of radius a; we write l/a = m. 

Suppose the rod placed in a uniform field H parallel to its axis, which 
we take as the axis of x ; then we have 

V = -Hr. (5) 

Our object is to calculate the functions Vi, V 2 ,... defined in (2). 



VOL. lxxix.—A. 


* Cf. lliecke, loc . cit ., p. 481. 


D 
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Using cylindrical co-ordinates (x, p, 0 ), we have in general 
_ f dV n (x f , p') da __ r 3V n da 


47 rY n+1 (x,p) = I ™ = I "15 w . . (6) 

J dv r ) dv >/(x'—x) 2 + p 2 + p' 2 — 2pp' cos0' 

Also the element of surface da is equal to p'ds d0\ where ds is equal to dx 
or dp according as the point (V, p) is on the curved side or a flat end of the 
rod; and putting 

TS? = f?+p'*-2pp' cos 0\ 


we have 


4ttY„ + i (x,p) = ['cM' jp' ds '-*>J 0 (XK) dX 


= 2tt Jp'-J ds (Xp) Jo (Xp') d\ 


( 7 ) 


where the integral with respect to ds extends over the path ABCD, and the 
upper or lower sign is taken in the exponential according as x*—x is 
negative or positive. 

Thus we find 

v n+1 (x, p) = * 9 dp' p-*(*+OJ 0 (Xp) Jo (V) 

-i p[ A ? p- x <'-*>Jo (V) J o (V) dX 

— 2 | a ^ X ' | « A(l '“^Jo (Xp) Jo (Xa) d\ 

— aAx ' | e -A ^' _:c) Jo (Xp) Jo (Xa) dX. (8) 

From symmetry in the case under consideration, this reduces to 

V .+1 (*,#>) = /(V p-"8inh(X®)Jo(X/>) Jo(V)^ 


f ( adx'[ e~ kx ' sirih (Xa:) J 0 (Xp) Jo (Xa) 
J x \ / P' = « Jo 


G?X 


— f acfa/ ( e -A:K cosh (X#') Jo (Xp) Jo (Xa) c?X. 

J 0 \ 3 p /?' = a J 0 


( 9 ) 


Thus at any stage in the evaluation of the successive potentials, V n+1 is 
expressed as the potential of a surface distribution due to V n ; the first 
integral in (9) gives the part due to the flat ends of the rod and the other 
two integrals that due to the curved sides. 

The potential Vi can be deduced from (9), or obtained directly as the 
potential of two uniformly charged circular discs situated at the ends of the 
rod. 
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Using the relation 


■ft 

p* Jo (f , p f ') dp f — 1 (xayjx, 


o 


we obtain 


.00 


Yi = H a I e~ u sinh (\&) J 0 (Xp) Ji (Xa) dX/X. 


( 10 ) 


o 


§ 4. Calculation of the Total Moment of the Rod. 


The magnetisation at any point of the rod is equal to — k (3 /dx, d/dy,d/dz) <£, 
where <£ is the total potential, and the moment of the rod can be calculated 
by considering the normal component of the polarisation over the surface 
of the rod. Writing M for the total moment, we have 


2iTKp(^\ dp — 2 f 27tkox(^\ dx, (11) 


M = —21 


fjXJx =>l J Q \dp /p = ft 


where the first integral gives the part due to the ends of the rod, and the 
second integral that due to the sides. 

Now cj> is given in terms of Y and the intermediary potentials Y n according 
to the series given in (3) and (4) ; thus from (11) we find that it is required 
to evaluate expressions of the form 



dx . 


( 12 ) 


ft 


We shall consider these expressions as far as U 2 , and obtain the important 
terms in their expansions in powers of 1/m, where m is the ratio of the 
length of the rod to its diameter, and is assumed to be large. 

Then since 


(13) 


Further, from the value of Yi given in (10) we have 



—47ra 2 H xdx e~ Kl sinh (Xx) J 1 2 (Xa) dX. (14) 


If in the second part we perform the integration with respect to x , and put 
p for Xa, we find that (14) reduces to 


Ui = 27ra 3 H f (V 2 ™^ HM-d/j,. (15) 



Now we have 
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Also, since m is large, the important terms of the second integral are 
obtained by substituting the value of the Bessel function for small values of 
the argument; thus as far as the term m~ 2 we have 


f dfjL=— 

I o Z*' 


8m 


(16), 


Hence to the same order we have 


U, = 2™»H(±-i 


= 2ira 2 m 


\?>Trm 8m 2 J 




Now we have to evaluate U 2 in a similar manner; we have from the- 
general expression given in (9), 

= ~\°fS)f' dp ' f” Xe_Ai cosh (^) J o (V) Jo (V)^ 

+ f (~J~ ) adx' [ \e~ u cosh N Xa/) J 0 (Xp) J 0 (Xa) dX 
J q\ Op /° J o 

= -*©,-* £©/«' -MV) 

+ | (^~rj a ^ x ' j ~Xe~ M cosh (Xx') J 0 (Xp) J 0 (Xa) dX. (18V 


And, similarly, from (9) we have 

av L 

■) = -M-5ri.+ i o 


I s ) ,r 


( 19 )* 


It may be noted that if the intermediary potentials Y n be regarded as due - 
to successive surface distributions cr n , then the equations (18) and (19) are 
examples of general relations between consecutive surface densities. 

Using the values of (18) and (19), we have 

U 2 = -l jVp (^i) ; dp-l 2irpdp^(~?jf'dp' |“xe- 2 « Jo (Xp) J 0 (Xp') dX 

+ 2Z | 2irpdp | adx'j cosh (Xx') J 0 (Xp) J 0 (Xa) dX 

— | 2ttcix dx + 2 | 2nraxdx J (^-^j^p'dp | Xe~ xl sinh (Xx) Ji (Xa) 

Jo (Xp ) dX. (20) 
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Using the formulae 

f d\ i\pcj) ( p ) J n (\p) Jn (Xr) = <£ (r), q < r < p, 

Jo J^ 

r 00 1 

I Jo (\/>) Ji (\&) d\ = p < a, 

Jo a 

we find that (20) reduces to 

U 2 = -jVaz l^jdx-^a f a (jtf)f' d P' e ~ U sinh OO J o (VO J 1 (*«) Y 

+ 47ra^| | e~ AJ cosh (\a;') Jo (Xa) Ji (\a) c£\. (21) 

We substitute now for Vi from (10) and write 

x = a %; p f =z arj; l = am ; ya = \a; pci = \'a . 

Hence we obtain 

rm r 00 

U 2 = 27rHa 3 sinh (yai£) Ji 2 (yai) dyai 

Jo Jo 

—47rHa 3 f TfdTj j* cosh (m/^i) Jo(ya?;) Ji(yai) eifyai 

Jo Jo 

f sinh (mya) J 0 (pv) J 1 (m) — 

Jo M 

—47rHa 2 Z f d%[ e~ m ^sinh (yai0 Ji 2 (yai)tZyai 

Jo Jo 

f e ~w/x j o j x dp ^22) 

Jo 


The first two terms give the following: 


dpu 

7 

dp, 

M ’ 


27rHa 3 w f e~ m ^ cosh (myai) Ji 2 (yai) —47rHa 3 f e~ m ^ sinh (mya) Ji 2 (ya) 

Jo Mi Jo 

/*1 <»GO *CO 

— 27rHa 3 J rjdr) j e~ 3m ^i j 0 (yai?;)J 1 (yai)^yai j sinh (mya) J 0 (ya? 7 ) Ji (ya) 

Jo Jo Jo 

Also, in the third term of (22) we perform the integration with respect 
to £; then in all the terms of U 2 we apply the previous method of approxi¬ 
mation, that is, we consider m large and substitute the lower terms of the 
expansions of the Bessel functions for small values of the argument. Using 
also the relation 

rww i„ = } , 

Jo M 

and keeping to the order in which the parts of U 2 are given above, the 
approximation indicated gives, to the order 1/m 3 , 
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U 2 = 7rHa 2 //l+ 1 


_8_ 1_1_1 J__1_\ 

2~*~16m 2 37 rm 4m 2 67m 3 2 87 rm 3 16m 2 / 


= ,rHa a J-JL + 1 


\ 87 rm 4?/i 2 247r?ft 3 


(23) 


Comparing U 2 with Ui given in (17), we see that they both begin with a 
term of the same order 1 /m; but we have 


Ui + U 2 = -i-irHa*/- 1 


24 


nr 


(24) 


Eeturning now to the series (3) and (4) for the induced potential Q, we 
see that in calculating the total moment of the rod in this manner it is of 
advantage to use the series 

Q = ^V 1 + te) 2 (Y 1 +Y 2 )+... ) fM = 1 + 4™. 

Then, if we denote by I m the mean moment per unit volume obtained by 
calculating the moment of the whole bar, we have from (13), (17), and (24)— 
as far as the order indicated— 



(25) 


§ 5. The Mean Induction through the Central Section. 


Another quantity which we require to calculate is the mean polarisation 
over the central transverse section of the bar; if this be denoted by I c , 
we have 


where 



<j> = — H& + Q. 


(26) 


To evaluate the contributions made by Vi and V 2 , we have from ( 10 ), 



= f pdpi e~ xl Jo (\p) Ji (Xtf) d\ 

a Jo Jo 

= 211 C e-^Ji 2 {fi)dfi/fi 

J 0 


= II 


Inr 



( 27 ) 


to the order indicated, using previous method of approximation. 
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From the general expression in (9) we have V 2 given in the form 

ra r® 

V 2 = —Ha p'dpx e _A ^ cosh (Xi/)J 0 (Xip')Ji (Xia) dXi 

Jo Jo 

f ee 

e~ xl sinh (X#) J 0 (Xp ) Jo (Xp') dX 

0 

-f Ha 2 j* dx'[ e~X l sinh (X\x') J 1 2 (Xia) dXi 

J z Jo 

f 0 _Az 'sinh(X#) J 0 (Xp) Jo(Xa)c£X 

J 0 

•f Ha 2 f rfa' f e~ x ' 1 sinh (Xi#') J 1 2 (Xia) dX 1 

J° J 0 

j* e~** cosh (Xx') J 0 (Xp) J 0 (Xa) dX. (28) 

J 0 


Hence we have 

av 2 \ 


-g— 2 ) = —Ha j p dp | e~^ 1 cosh (Xd) J 0 (Xip') Ji (Xia) dXi 

f Xe~ xl J 0 (X/>) J 0 (X//) dX 

Jo 

+ Ha 2 f dx' j* e~X l sinh (X^') J 1 2 (Xia) dX 1 

J o J° 

Xe~ Xx> J 0 (Xp) J 0 (Xa) dX. (29) 
. 0 

The mean value of this taken over the circle x = 0, p = a, is given by 

= —2H f p dp f c _A ^cosh (X x l) J 0 (Xip') Ji (Xia) dX x f e~ A *Ji(Xa) J 0 (Xp')dX 

Jo Jo 


+ 2Ha 


n r® 
dx 6 _A ^si 
Jo Jo 


sinh (Xi^')Ji 2 (Xia) d\\ 


e Xx> Jo (Xa) Ji (Xa) dX. (30) 


We proceed to evaluate these integrals for m large by substituting the 
expansions of the Bessel functions as before; then the first integral in (30) 
becomes, to the required order, 


= _tj / _1 _13 _ 1 _\ 

\4m 2 32 to 4 / ' 


(31) 


The second integral in (30) is equivalent to 

e -mfx e -'>n(tx + ‘2fj. 1 ) g-wiMi 


H 


ca 


- + 


/ii —ya p\-\-p P'l — P' A 6 ! J 


Mi 2 (/ii) Jo(/i)Ji (f) dp, dp,^ 
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Substituting values of the Bessel functions of /*, or fix, or of both, according 
to the exponential term, this reduces to 

H (s=rS=>)- < 32 > 

Adding this to the part given in (31), we see that the mean value of 
dV 2 /dx over the central section begins with a term —23H/32m 4 , while the 
mean value of dYi/dx begins with a term of the order 1/m 2 . 

Then, substituting in (26) the values calculated above, we have, as far as 
the order shown, 

I. = J- 4 )-( W H jy 

")=>}■ < 33 > 

§ 6. Application to Dielectric Polarisation. 

If an ellipsoid of specific inductive capacity K is placed in a uniform 
electric field, it is known that the field within the ellipsoid is uniform, and 
the polarisation at any point is given by a formula 

I = k (H—NI); k = (K-1)/4t r. (34) 

The factor N is a numerical factor independent of I, and can be calculated 
theoretically; for a prolate spheroid, of which m is the ratio of the longer 
to the shorter axis, we have 

N = J=i{^3=! ll>s * ( “ +v/s=I) - 1 }' < 36 > 

the axis of revolution being along the direction of the field. Also, for an 
infinite rod of circular section, N is zero for a uniform field parallel to the 
axis of the cylinder, and is equal to 2 7r for a field at right angles to 
the axis. 

In the case of a rod of finite length placed in a uniform field, there is no 
longer uniformity of the field within the rod, but as we suppose the specific 
inductive capacity to be independent of the field, we have still an equation 
of the form (34). Bor a field parallel to the axis of the rod, we have from 
(25) to the first order of 1/m, the ratio of diameter to length, 

N = 16/3Km. (36) 

For a field at right angles to the axis we may still take N to be equal to 
2 t r. Now, in determinations of the specific inductive capacity of a dielectric 
in the form of a cylindrical rod, it is necessary to measure the couple 
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required to keep the rod at a given angle <jf> to the direction of the external 
field. If we put Y for the volume of the rod and use the values of N 
obtained above, we have 

Couple on rod = k ?Y ■ ^Tr—lG/SKm jp j s j n ^ cos ^ (37) 

r (l + 27r/c)(l + 16/c/3Km) r ■ v 

where k = (K—1)/47t. 


§ 7. The Demagnetising Factor of a Finite Bod , 

When the problem is one of magnetic induction, it is clear that the 
previous results only supply an approximation to a solution; for, since the 
field is not uniform within the rod, it follows that the magnetic permeability 
jjb varies along the rod, thus complicating the problem to be solved. In fact, 
the demagnetising force will not be expressible in the form NT, where N is 
-a numerical factor independent of I. 

On the supposition that there is such a constant N, at least approximately, 
for long thin rods, Du Bois* * * § made a comparison between the values of N 
for ovoids and cylindrical rods of the same dimension-ratio m ; the factors 
for the rods were calculated from experiments by Ewing for the longer rods 
.and by Tanakadate for the shorter ones. 

It was pointed out by Mannf that in the former experiments ballistic 
methods were used, but in the latter magnetometric methods. Consequently 
the quantities measured were those denoted by I m and I c above; and the 
factor N deduced by the two methods will, in general, be different, and may 
be denoted by N w and N$, being the magnetometric and ballistic values of 
the demagnetisation factor. From the series given in (25) and (33) we have, 
for large values of m, the approximate values 


N w 


16 . 
3 gm 5 




2ir 

m 2 * 


(38) 


Benedicks^ further investigated the values of N w and N& for rods and, 
in particular, their variation with the intensity of magnetisation; the value 
of the demagnetising force N$I at the centre of the rod was found to increase 
with the field, and then after a certain point to decrease rapidly. 

Using the ballistic method, Lamb§ determined the induction through 
sections taken at intervals along a bar, and was thus able to trace curves 


* Du Bois, 4 Annalen der Physik/ vol. 46, p. 495, 1892. 

+ C. B. Mann, ‘ Physical Beview,’ vol. 3, p. 359, 1895. 

f Benedicks, 4 Annalen der Physik,’ vol. 6, p. 726, 1901. 

§ C. G. Lamb, ‘ Phil. Mag./ vol. 48, p. 262, 1899. 
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showing the distribution of induction and force along the bar and the 
consequent effective distribution of magnetic permeability. 

Searle and Bedford* investigated the demagnetising force h at the centre 
of a long cylindrical rod. Determining the values of Ji for increasing values 
of the external force H, they found that h, increases at first more rapidly 
than 27 rl c /m 2 , but after passing a maximum it descends to a minimum 
somewhat less than that value; the latter result is interpreted as being due 
to hysteresis effects. 

A quantity measured in several experiments was the ratio I OT /I C , giving 
the effective length of the bar. 

Using expressions (25) and (33) above, we have approximately 



With a rod for which m was 25, Benedicks found this ratio increased to 
0*84 at the highest value of H used in the experiment. 

Reviewing the experimental results, it appears that in the case of magnetic 
induction the facts of the non-uniformity of the force within the rod and 
the consequent variation of effective susceptibility are too important to be 
neglected; hence, in this case, the calculations of the previous sections 
supply limiting values towards which the quantities approach as the 
magnetisation tends to saturation. 


* Searle and Bedford, ‘Phil. Trans., 5 A, vol. 198, p. 98, 1902. 










